New hadrons as ultra-high energy cosmic rays 



m 
o 
o 

(N 
£ 

(N 

(N 
> 

VO 

(N 
O 

m 
O , 

^ : 

i 

<D . 
43 ■ 



X 



Michael KachelrieB a , Dmitry V. Semikoz a ' & and Maria A. T6rtola c 
1 Max- Planck- Institut fur Physik (Werner-Heisenberg-Institut) , Fdhringer Ring 6, D-80805 Miinchen, Germany 
b Institute for Nuclear Research of the Academy of Sciences of Russia, Moscow, 117312, Russia 
c Instituto de Fisica Corpuscular - C.S.I.C., Universitat de Valencia, 
Edificto Institutos de Paterna, Apt. 22085, E-46071 Valencia, Spain 
(Dated: February 18, 2003) 

Ultra-high energy cosmic ray (UHECR) protons produced by uniformly distributed astrophysical 
sources contradict the energy spectrum measured by both the AGASA and HiRes experiments, as- 
suming the small scale clustering of UHECR observed by AGASA is caused by point-like sources. In 
that case, the small number of sources leads to a sharp exponential cutoff at the energy E < 10 20 eV 
in the UHECR spectrum. New hadrons with mass 1.5-3 GeV can solve this cutoff problem. For 
the first time we discuss the production of such hadrons in proton collisions with infrared/optical 
photons in astrophysical sources. This production mechanism, in contrast to proton-proton colli- 
sions, requires the acceleration of protons only to energies E < 10 21 eV. The diffuse gamma-ray 
and neutrino fluxes in this model obey all existing experimental limits. We predict large UHE neu- 
trino fluxes well above the sensitivity of the next generation of high-energy neutrino experiments. 
As an example we study hadrons containing a light bottom squark. This model can be tested by 
accelerator experiments, UHECR observatories and neutrino telescopes. 

PACS numbers: 14.80.-j, 98.70.Sa 



I. INTRODUCTION 

UHECRs with energies above 10 19 eV have been ob- 
served in all relevant experiments, i.e. Volcano Ranch fil, 
Haverah Park 0, Fly's Eye 0, Yakutsk ||, AGASA Q 
and HiRes Their arrival directions are distributed 
uniformly over the sky without significant correlation 
with the galactic or the super-galactic plane. This 
isotropic distribution is consistent with the simplest 
model for UHE primaries, in which protons are accel- 
erated in extragalactic, uniformly distributed astrophys- 
ical sources. However, UHE protons with energies above 
E > 4 x 10 19 eV interact with cosmic microwave back- 
ground (CMB) photons and lose quickly energy through 
pion production within 50 Mpc. As a consequence, a cut- 
off in the UHECR spectrum, predicted already in 1966 
by Greisen, Zatsepin and Kuzmin (GZK) pfl, should show 
up for uniformly distributed sources at ~ 5 x 10 19 eV. 
This cutoff is not observed by the ground array experi- 
ment with the largest exposure, AGASA, while the first 
monocular results of the HiRes fluorescence telescope are 
in agreement with the GZK cutoff. The exposure at the 
highest energies of all other experiments is too small to 
allow for a definite conclusion about the presence or ab- 
sence of the GZK cutoff. 

Fortunately, the next generation experiment Pierre 
Auger Observatory 8] which is a combination of an array 
of charged particles detectors with several fluorescence 
telescopes is currently under construction. It will not 
only be able to resolve possible systematic differences be- 
tween the ground array and fluorescence telescope tech- 
niques, but will also increase the statistics of UHECR 
data by an order of magnitude. The telescope array 
project, also based on the fluorescence technique, may 
serve as the optical component of the planned northern 



Pierre Auger site 9] . There are also plans for sp ace based 
observatories such as EUSO [13 and OWL [0 with even 
bigger acceptance. 

Assuming that the GZK cutoff will be confirmed by 
future experiments does not resolve the UHECR puzzle. 
Since all experiments including HiRes see events with 
energies E > 10 20 eV, their sources should be located 
within the distance R < 50 Mpc. Otherwise, the GZK 
cutoff is extremely sharp and in contradiction even to 
the UHECR spectrum measured by HiRes (cf. Fig. |BJ . 
But there are not many astrophysical sources within this 
distance from the Earth known able to accelerate par- 
ticles to the highest energies. Moreover, these sources 
are not located in the directions of observed events. An- 
other problem is a statistically significant (4.6 a for en- 
ergies above 4 x 10 19 eV) clustered c omp onent in the 
arrival directions of AGASA data [H El Q [l5|. The 
sensitivity of the other experiments for clustering at the 
energies E > 4 x 10 19 eV is much smaller, either be- 
cause of the smaller exposure at the highest energies 
(Yakutsk) or because of a poor two-dimensional angu- 
lar resolution (HiRes in monocular mode). At lower 
energies 10 19 eV < E < 4 x 10 19 eV a clustered com- 
ponent still exists in the AGASA data but with a 
reduced significance of 2.3 a. The Yakutsk experiment 
also observes a clustered component in the energy region 
£>2.3x 10 19 eV with a chance probability 2 x 10~ 3 or 
~ 3cr using Gaussian statistics |13|. 

The puzzle of the GZK cutoff can be solved in two 
different ways. The first one supposes that the sources 
of UHECR are located nearby. Then the extragalactic 
magnetic field should be strong enough, B > 0.3^G, to 
deflect UHECR with E > 10 20 eV and magnetic lensing 
could be responsible for the clustered component ^3 ■ A 
problem of this solution is the difficulty to construct real- 
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istic maps of the matter and magnetic field distribution 
in the nearby Universe. Simulations done so far repro- 
duced the energy spectrum and the clustered component 
assuming 10-100 sources, but without using realistic lo- 
cations of the sources. Another difficulty is that magnetic 
lcnsing, although reproducing the clustered component, 
predicts in general a broad angular distribution of this 
component, while the data are within the experimental 
angular resolution. This solution is even more problem- 
atic if future experiments show the absence of the GZK 
cutoff: fixing the luminosity of the local sources to the 
UHECR flux above the GZK cutoff then results in a too 
large UHECR flux at lower energies, where sources from 
the whole Hubble volume contribute. 

The second way is to suppose that the clustered com- 
ponent is due to a neutral particle which is not deflected 
by (extra-) galactic fields. In this case one can look for 
correlations of UHECR arrival directions with astrophys- 
ical objects. Tinyakov and Tkachev have recently found a 
significant (more than 4 a) correlation with BL Lacs • 
The BL Lacs which correlate with the UHECRs are lo- 
cated at very large (redshift z ~ 0.1) or unknown dis- 
tances. If it can be shown with an increased data set 
of UHECRs that this correlations holds also at energies 
E > (6—10) x 10 19 eV, then protons alone can not explain 
the UHECR data and a new component in the UHECR 
spectrum is needed. 

The simplest possibility is that this new component 
is due to extremely high energy (E > 10 23 eV) photons 
emitted by distant sources. They can propagate several 
hundred Mpc constantly losing energy and thereby creat- 
ing secondary photons also inside the GZK volume [l8| . 
However, this model requires extremely small extragalac- 
tic magnetic fields, B < 10~ 12 G, and the minimal pos- 
sible radio background. Besides, one needs to accelerate 
protons to E > 10 24 eV in order to create such photons. 
An acceleration mechanism to these extreme energies is 
not known. 

Another possibility is that the events beyond the GZK 
cutoff are related to the Z burst model [l^. In this 
model, UHE neutrinos interact with the relic neutrino 
background producing via the Z resonance secondary 
protons and photons. The big drawback of this scenario 
is the need of an enormous flux of primary neutrinos that 
cannot be produced by astrophysical acceleration sources 
without overproducing the GeV photon background [2jj . 
Also in this model, primary protons have to be acceler- 
ated to extremely high energies, E > 10 23 eV, in order 
to produce E = 10 22 eV neutrinos. 

Conventionally, acceleration mechanisms allow to ac- 
celerate protons in astrophysical sources only up to E < 
10 21 eV. If one considers this maximal energy as a se- 
rious upper limit, both possibilities discussed above are 
excluded and some kind of new particle physics beyond 
the standard model is required. The most radical op- 
tion is violation of Lorentz invariance |2l| . A more con- 
servative, though for same tastes still too speculative, 
possibility are decaying super-heavy relics from the early 



Universe |22| . This model cannot explain the correlation 
of UHECR arrival directions with BL Lacs, and could 
be excluded if these correlations are found also at en- 
ergies E > (6 — 10) x 10 19 eV. In theories with a fun- 
damental scale of gravity as low as O(TeV), the inter- 
actions of neutrinos with nucleons are enhanced com- 
pared to the standard model by the exchange of Kaluza- 
Klein gravitons [2j| or the production of black holes [24| . 
In Refs. [2^1 it was claimed that neutrinos could have 
hadron-like cross sections at UHE and be responsible for 
the observed vertical air showers. However, unitarization 
slows down the growth of the neutrino-nucleon cross sec- 
tion [2|| and moreover, in the case of exchange of 
Kaluza-Klein gravitons, the energy transfer to the tar- 
get nucleon is too small ; 26J. Thus the neutrino-nucleon 
cross-section can be up to a factor 100 larger than in the 
standard model, but neutrinos in both cases still resem- 
ble deeply penetrating particles and cannot imitate air 
showers initiated by nucleons. 

Another possibility is that new particles are directly 
produced in astrophysical sources. A model with an 
axion-like particle, i.e. a scalar which can mix with a 
photon in the presence of external magnetic fields, was 
suggested in Ref. [2^]. Axion-like particles can be also 
produced by photons emitted by astrophysical sources 
via axion-photon oscillations [29| . In supersymmetric 
(SUSY) theories with conservation of R parity, the light- 
est supersymmetric particle (LSP) is stable and can be 
a HE primary. This possibility was for the first time 
seriously discussed in connection with Cyg X-3 in the 
80's [30L l31j. More recently, the production and inter- 
actions of both the neutralino and the gluino as LSP at 
UHE were examined in Ref. j33| ■ The authors concluded 
that only a light gluino could be produced in reasonable 
amounts by astrophysical accelerators. Reference [33] 
calculated the neutralino production in proton-proton 
collisions and found that neutralinos produced in astro- 
physical sources cannot be an important UHE primary: 
Since the production cross section of neutralinos is too 
small, this model predicts either a negligible flux of UHE 
neutralinos or is not consistent with measurements of the 
diffuse gamma-ray background [34| and with existing lim- 
its on the neutrino flux at ultra-high e nerg ies. The latter 
limit were obtained by the Fly's Eye Hj, AGASA [36| . 
RICE H3, and GLUE [13 experiments. 

SUSY models with a strongly interacting particle as 
LSP or next-to-lightest SUSY particle (NLSP) are more 
interesting for UHECR physics. Hadrons containing a 
gluino were first suggested by Farrar as UHECR pri- 
mary [39I l4(i| . Her model, a light gluino g together with 
a light photino such that the photino could serve as cold 
dark matter candidate, is excluded 0,^2- Motivated 
by the correlation of UHECR with BL Lacs, the produc- 
tion of light gg bound-states in astrophysical accelerators 
was suggested in Ref. However, the light gluino win- 
dow seems to be now closed by Ref. also for generic 
models with a light gluino. 

In this paper, we start from a model- independent, 
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purely phenomenological point of view. Since the ob- 
served extensive air showers (EAS) are consistent with 
simulated EAS initiated by protons, any new primary 
proposed to solve the GZK puzzle has to produce EAS 
similar to those of protons. An experimentally still open 
possibility are photons as UHECR primaries: at 90% 
C.L., - 30% of the UHECR above E > 10 19 eV can 
be photons ^5|- However, the simplest possibility consis- 
tent with air shower observations is to require that a new 
primary is strongly interacting. The requirements of ef- 
ficient production in astrophysical accelerators as well as 
proton-like EAS in the atmosphere ask for a light hadron, 
< 3 GeV, while shifting the GZK cutoff to higher ener- 
gies results in a lower bound for its mass, > 1.5 GeV [43j |. 
From these requirements, we derive general conditions on 
the interactions of new UHE primaries. 

As specific example, we investigate the case of a bot- 
tom squark containing hadron which we call "shadron" 
from now on. Our conclusions are however indepen- 
dent of the underlying particle physics model for the 
shadron. The required properties of the shadron will be 
parameterised as function of its mass, production and 
interaction cross-sections. Other suitable realizations of 
shadrons could be new (meta-) stable hadronic states like 
H-dibaryons ^(| . Another candidate could be connected 
to the exotic, charged hadronic state with mass 2.3 GeV 
discovered recently by BaBar which was suggested 
by the experiment as a four-quark state. If this sugges- 
tion would be true, then this four-quark state could be 
related to a new metastable neutral hadron. 

We find that proton-proton collisions in astrophysical 
accelerator cannot produce high enough fluxes of new pri- 
maries without contradicting existing measurements of 
photon [34j and neutrino fluxes j3j| |3(| 13]} By contrast, 
we find for a light shadron with mass < 3 GeV and the as- 
trophysically more realistic case of UHE proton collisions 
on optical/infrared background photons no contradiction 
with existing limits. Also, the required initial proton en- 
ergy is not too extreme, E < 10 21 eV, which is com- 
patible with existing acceleration mechanisms. The only 
essential condition for the sources is that they should be 
optically thick for protons in order to produce these new 
hadrons. (This condition is similar for all models with 
new particles produced by protons) . Below we will show 
that at least some of BL Lacs correlated with UHECR 
obey this condition. 

One of the important features of the proposed model, 
and any model in which the production cross section 
tjpy^s °f a n ew particle S is much smaller than the to- 
tal proton-photon cross section cr p7 , is the high flux of 
secondary high-energy neutrinos. This neutrino flux is 
connected via the relation Fq^u p1 / a pl —,s to the max- 
imal contribution of S particles to the cosmic ray flux, 
Fqk ~ I/-E20 eV/(cm 2 s sr). It can be detected by future 
UHECR experiments like the Pierre Auger Observa tory 
@, the Telescope Array @, EUSO [13 and OWL [nj. 
Alternatively, such neutrino fluxes can be detected by 
triggering onto the radio pulses from neutrino-induced 



air showers |48j . Acoustic detection of neutrino induced 
interactions is also being considered ^3 ■ There are plans 
to construct telescopes to detect fluorescence/ Cerenkov 
light from near- horizontal showers produced in mountain 
targets by neutrinos at intermediate energies 
Moreover, if the sources are optically thick for protons, 
the neutrino flux can be significant both at high en- 
ergies and down to energies 

lQie-17 eVj depending on 
the pion-production threshold on optical/infrared pho- 
tons [52). Therefore, one may observe neutrinos from the 
same sources both by future UHECR experiments and by 
neutrino telescopes like AMANDA [HI, ICECUBE [Hi. 
GVD 55], ANTARES 0, NESTOR~|57| or NEMO 0. 

The paper is organized as follows. We start with a 
discussion of the spectrum of UHECR protons produced 
by a small number of extragalactic astrophysical sources 
in Sec. [H] Then we consider models containing light 
strongly interacting particles, shadrons, and their sta- 
tus. In Sec. EH we discuss the propagation of shadrons 
through the Universe. Their interactions in the atmo- 
sphere were investigated in detail before, so we shall just 
briefly recall the main characteristics in Sec. El Sec- 
tion El is devoted to a detailed analysis of shadron pro- 
duction in astrophysical sources. In Sec. IVIll we discuss 
all astrophysical constraints which shadrons have to obey 
to be viable UHECR primaries. In Sec. IVIIII we discuss 
the particular case of BL Lacs as sources of UHECRs. 
Finally, we summarize our results in Section llXI 



II. PROTONS FROM UNIFORMLY 
DISTRIBUTED SOURCES 

The HiRes experiment published recently their data 
from monocular observations They showed that the 
UHECR flux is consistent with the GZK cutoff expected 
for uniformly, continuously distributed sources. As a 
result, the simplest model of UHECR — protons accel- 
erated in uniformly distributed, extragalactic sources — 
seems to be a convincing explanation of their data. The 
authors of Ref . [H^ found as fingerprints of the expected 
interactions of UHE protons with CMB photons a dip 
at E ~ 1 x 10 19 eV, a bump [§3] and the beginning 
of the cutoff in the measured spectra of four UHECR 
experiments. The agreement of the spectral shape cal- 
culated for protons with the measured spectra is excel- 
lent, apart from an excess in the AG AS A data above 
i?>8x 10 19 eV. These findings point to an AGN ori- 
gin of UHECR below E < 10 20 eV and to protons as 
primaries. Despite the fact that the AGASA experiment 
sees a significant number of events above the GZK cut- 
off 5] , the model of proton primaries from extragalactic 
sources looks very attractive, because it does not require 
new physics. 

The model of uniformly, continuously distributed 
sources is based on the assumption that the number of 
UHECR sources is so large that a significant fraction of 
sources is inside the GZK volume. However, as it was 
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shown in a number of works 0], j^], the small 

scale clustering of UHECR observed by AGASA allows 
to estimate the number of UHECR sources assuming that 
their distribution and luminosity is known. For the sim- 
plest model of uniformly distributed, similar sources their 
number is about several hundreds, 200 - 400. If we dis- 
tribute these sources uniformly in the Universe, the num- 
ber of sources in the GZK volume with R = 50 Mpc 
would be of the order of 10~ 3 . This would mean that the 
nearest source should be at the redshift z = 0.1. 

A more conservative and self-consistent estimate uses 
the fact that protons with energies E > 4 x 10 19 eV 
observed on Earth can pr opag ate at most from redshift 
z = 0.2 (see, e.g., Fig. 2 in (13 )• Distributing the sources 
within a sphere at z = 0.2 around the Earth, the closest 
source is at the distance R = 100 Mpc. Note also that in 
the particular case of BL Lacs as UHECR sources, which 
we discuss in Section IV 1111 the closest potential sources 
are at redshift z ~ 0.03. 

We show now that the statement that UHECRs with 
E > 10 20 eV are protons from nearby sources is in contra- 
diction to the total number of sources estimated including 
events below the GZK cutoff. Using Poisson statistics, 
the total number of sources S is fixed by the number of 
observed singlets N\ and doublets N 2 , 



Ni 
N 2 



Sn , 

S^ 
b 2 



(1) 
(2) 



where n is the average number of events from a given 
source and we assumed equal flux from all sources. From 
Eqs. ||TJ and J5J one obtains the number of sources 
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As it was shown in Ref. [6IJ , the value Eq. © is a model 
independent lower bound on the number of sources for 
given values of Ni and N 2 . In the model of homoge- 
neous distribution of sources with equal luminosity, the 
estimate for the number of sources becomes roll. 
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where N to t is the total number of observed events and 
iV c i = 2N2 is the number of events in clusters. In Eq. 
it is assumed that N to t 3> N c \. Note that therefore 
Eq. Q gives always a larger number of sources than 
Eq. J3Jl because of the extra factor N% Q t/N c i 3> 1. 

We estimate next the number of sources assuming that 
all UHECR with £>4x 10 19 eV are protons. Follow- 
ing [H i2, we use 14 events with E > 10 20 eV and 
one doublet. Calculating S with Eq. (0) gives S ~ 100 
as a minimal number of sources, while the more realistic 
Eq. J2J gives S ~ 700. If we apply the same analysis to 
the AGASA data [3 with E > 4 x 10 19 eV, we have 



N2 = 6 (we count for simplicity the triplet as one dou- 
blet) and Fi = 46. Then Eq. © gives S ~ 176, while 
Eq. Q gives S ~ 1200. 

Are these two estimates consistent with the idea that 
all UHECR are protons? To answer this question, we 
calculate the expected number of proton sources using 
E > 4 x 10 19 eV when there are S ~ 100 — 700 sources in 
the GZK volume. Protons with E ~ 4x 10 19 eV can reach 
us from z = 0.2, or R tot ~ 1000 Mpc. Conservatively 
assuming that all events with E > 10 20 eV come from 
within the GZK distance R = 50 Mpc fin |6lll62T| R = 25 
Mpc was used), we obtain with Eq. © as expected num- 
ber of sources S to t = 0Rtot/-RGZK) 3 x 100 = 8 x 10 5 . 
Using instead Eq. Q , the expected number of sources is 
Stot = 5.6 x 10 6 . These estimates should be compared to 
to the ones from AGASA clustering data, Sagas a ~ 176 
or Sagas A ~ 1200, using Eq. or Eq. (QJ respectively. 
Since the Poisson probability to observe Sagasa instead 
of 

Stot events is practically zero, the chance probabil- 
ity to obtain these two event numbers is equal to the 
chance probability of clustering. We conclude therefore 
that the model in which all UHECR with E > 4x 10 19 eV 
are protons from uniformly distributed point sources is 
inconsistent with the small scale clustering observed by 
AGASA. 

One can argue that 14 UHECR events with E > 
10 20 eV is an optimistically high number and that the real 
number of such events is much smaller because the exper- 
iments estimate wrongly the energy of UHECR events. 

We conservatively take only the four highest energy 
events from all experiments, including one Fly's Eye 
event, two AGASA events and one HiRes event. In this 
case we have 4 single events and no doublets. We can 
estimate the number of sources from absolute minimal 
bound Eq. © if we assume that the average number of 
doublets is less then one, i.e. e.g. N2 — 0.5. Then there 
are S = 16 sources in the GZK volume with -Rgzk = 50 
Mpc. Again, in a volume with i? to t ~ 1000 Mpc there 
are S ~ 128.000 sources, in comparison with up to 1200 
required by AGASA data above E >4x 10 19 eV. 

Thus, if the clustered component in the AGASA events 
with energy E > 4 x 10 19 eV is due to point-like sources, 
the expected number of sources is of the order of several 
hundreds up to 5 ~ 1200, depending on the estimate 
used. These sources are distributed in a volume with 
Rtot ~ 1000 Mpc. Assuming that the UHECR events 
with E > 10 20 eV are protons requires 10 — 400 sources 
in the GZK volume with -Rgzk = 50 Mpc. These two 
facts are in contradiction, if all UHECR are protons. In 
other words, if UHECR with E > 4x 10 19 eV are protons, 
we should have less than one source, namely S < 0.1, in 
the GZK volume. 

Let us now discuss the consequences of a small number 
of sources for the model of uniformly, continuously dis- 
tributed point sources of protons. For our calculations, 
we have used the code developed in Ref. [6^ , in which all 
important effects (pion production, e + e~ production and 
the expansion of the Universe) are taken into account. 
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FIG. 1: UHECR flux measured by the HiRes experiment |(J. 
The thin solid line corresponds to an uniform, continuous dis- 
tribution of proton sources in the Universe with emission spec- 
trum 1/E 2 - 7 and _B max = 10 21 eV. The dotted curve is for the 
same model, but with no sources within 50 Mpc from the 
Earth. The thick solid line corresponds to no sources within 
Zmin = 0.03, the dashed line to z m i n = 0.1. 
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FIG. 2: UHECR flux measured by the AGASA experi- 
ment [f| . All other parameters are the same as in Fig. 



Essentially, we have repeated for the case of proton pri- 
maries and BL Lacs as sources the calculations made 
in 18] for photons in more details. In Figs. [I] and we 
show with thin solid lines the spectrum of continuously 
distributed sources of protons with emission spectrum 
1/E 2 - 7 and £ max = 10 21 eV as in Ref. |Hj|. The dot- 
ted, thick solid and dashed lines are for the same model, 
but with no sources within 50 Mpc, z m j n = 0.03 and 
•Zmin = 0.1 around the Earth, respectively. The mini- 
mum distance of z m - ln — 0.03 corresponds to the BL Lac 
distribution. 

Let us concentrate on Fig.^ which shows the measured 
spectrum of HiRes and where the fit model of [EjJ with 
an infinite number of sources (thin solid line) works well. 



If there are no sources within 50 Mpc (dotted curve), the 
two highest HiRes data points are well above the model 
fit. For the BL Lac case where the closest known sources 
are at z m i n = 0.03 two additional experimental points are 
away from the fit. Finally, for an uniform distribution 
of 400 proton sources over the Universe, or z m ; n = 0.1, 
the disagreement above the cutoff becomes even worse. 
Note that we are only concerned about energies above 
~6x 10 19 eV; at lower energies, the quality of the fitted 
model can be easily improved by a readjustment of the 
fit parameters. The same figure with experimental data 
from AGASA is shown in Fig. El 

Thus, if the clustered component of the AGASA data 
for E > 4 x 10 19 eV (which has a statistical significance 
of 4.6cr) is not a statistical fluctuation or the result of 
magnetic lensing, the expected relative small number of 
UHECR sources is inconsistent with the model of proton 
primaries emitted by uniformly continuously distributed 
sources both for the HiRes and AGASA data. This 
means that both the AGASA and HiRes data require the 
introduction of a new component (not protons) in the 
UHECR spectrum. In the following sections we will con- 
sider new light hadrons with the mass of 2-3 GeV as such 
a new component. 



III. LIGHT STRONGLY INTERACTING 
PARTICLES: MODELS AND EXPERIMENTAL 
STATUS 

We shall concentrate our discussion on the definite case 
of the light sbottom, which is rather predictive and al- 
lows various tests of our assumptions. Since an useful 
UHE primary should be stable or quasi-stable with life- 
time t > 1 month (cf. Sec. II V|) . only the LSP or the 
NLSP are possible candidates as new UHE messengers 
in SUSY models with (approximately) conserved R par- 
ity. The NLSP as UHE primary can be realized if it 
has a very small mass splitting with the LSP or if the 
LSP is the gravitino; in the latter case the gluino decays 
gravitationally and its lifetime can be long enough. 

Theoretically the best motivated candidates for the 
LSP are the neutralino x an d the gravitino G3/2. While 
in minimal supergravity models the LSP is the light- 
est neutralino (in some part of the parameter space it 
is the sneutrino), in models with gauge- mediated SUSY 
the LSP is normally the gravitino. Recently, a light bot- 
tom squark b with mass ~ 2 — 6 GeV has been sug- 
gested 0, motivated by the large bottom quark pro- 
duction cross section measured at the Tevatron [(55| : The 
long-standing puzzle of overproduction of bb pairs can be 
solved if there exists additionally to a light sbottom a 
light gluino with mass m~ g ~ 12 — 16 GeV. 

Bottom squarks as LSP can either form charged B~ = 
(bu) or neutral B° = (bd) (plus charged conjugated) two- 
quark states. Since qu states are generally lighter than qd 
states, it is likely that the charged = (bu) is lighter 
than the neutral B° — (bd). But their mass difference 
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will be very small, e.g., for the usual B system m B a — 
m B ± ~ 0.33 ± 0.28 MeV, and we consider therefore the 
question if the lightest state would be charged or neutral 
as open. Moreover, the mass difference in the B system 
could be smaller than the electron mass, and weak decays 
therefore kinematically forbidden. In this case, both the 
B^ and the B° would be stable. Apart from the two- 
quark states, there will be again baryonic three-quark 
states, like e.g. bud. These baryonic states can decay 
into a baryon and a B if kinematically possible. 

Theoretically, the lightest hadronic state could be elec- 
trically charged. Is it possible that a light, stable charged 
hadron evaded detection? At Serpukov and the CERN 
ISR several searches for such particles were performed in 
the 70's [7l|7llzil- For example, the CHLM experiment 
excluded the range m > 2.4 GcV for stable hadrons with 
charge q = 1 [75|. Below 2.4 GeV, the production of an- 
tideuteron could hide other hadrons with a similar mass. 
Since the ratio R of antideuteron to pion production in 
these experiments is rather high, R ~ 5 x 10 -4 [7q|. and 
the mass resolution of these experiments not too fine, 
a significant fraction of deuterons could be misidentified 
stable charged hadrons. Also the TRISTAN experiment 
did not include the deuteron region in their search for 
massive stable hadrons [Tlj • While the LEP experiments, 
in particular DELPHI, could exclude generally charged 
shadrons down to masses of 2 GeV, the limit is in the 
case of a sbottom with small couplings to the Z boson 
weakened and sbottoms with masses below 5 GeV are al- 
lowed |73 ■ The ALEPH exclusion limit was not extended 
to masses below 5 GeV |79j . The CLEO experiment was 
able to exclude charged hadrons with mass m < 3.5 GeV, 
however only for fractionally charges 80]. Since there is 
however also no positive evidence for a stable charged 
hadron, we consider mainly the option that the lightest 
sbottom containing hadron is neutral. 

Next we recall that a light sbottom quark is consistent 
with electroweak precision observable and with the LEP 
higgs mass limit |8lj. Reference showed that this 
scenario implies a light stop, mj < 98 GeV, offering the 
TEVATRON run-II experiments the possibility to (dis-) 
prove indirectly the light sbottom case. The observa- 
tion of a bb resonance in e + e~ annihilation is difficult to 
extract from the background because the bb resonance 
has to be produced in a p-wave [83 ■ Its contribution to 
e + e~ — > hadrons is small compared to the error of these 
measurements. Since the lightest B behave as a stable 
particle in any accelerator experiment, their identifica- 
tion would require a dedicated analysis. We consider 
therefore a sbottom with mass 1.5-3 GeV as a viable op- 
tion and shall investigate its use as a UHE primary in 
the subsequent sections. 

Finally, we discuss the case of rather short-lived 
shadrons. Possible decays are b — > b + G 3 / 2 in models 

where the gravitino is the LSP or decays like G — ► tt + v, 
etc., if i?-parity is violated. In Ref. 66], it was argued 
that these decays can be excluded by proton decay exper- 
iments. If the life-time of g or b is close to the required 



lower limit of ~ 1 month, the shadrons produced by cos- 
mic rays and contained in the detector material have time 
to decay during the build-up and start phase of the exper- 
iment. Since these experiments are deep underground, 
they are well shielded and cosmic rays or shadrons can- 
not reach the detectors. In the case of detector using 
purified detector material, originally contained shadrons 
could be also extracted in the purification process, de- 
pending on the chemical properties of the shadrons. 

Light sbottoms do not contradict cosmological limits: 
The relative abundance of gluinos is ng/n 7 ~ 10~ 20 — 
10 _17 (mg/GeV) [72 and possible decays do not disturb 
Big Bang Nucleosynthesis (BBN) 0] • If baryon number 
also resides in baryons N containing sbottoms, then their 
number is suppressed by n^/riB ~ exp(— M^/Tqcd) ~ 
exp(— 3/0.16) ~ 10~ 7 , where Tqcd is the temperature 
of the QCD phase transition. Again there is no conflict 
with BBN. 

We should stress that light sbottoms in SUSY theories 
serve as merely as examples. Any particle physics model 
which has a (quasi-) stable particle with mass 1.5-3 GeV 
and interacts strongly with protons should have similar 
consequences for the physics of UHECR. 



IV. PROPAGATION THROUGH THE 
UNIVERSE: HOW TO AVOID GZK CUTOFF 

Strongly interacting particles S propagating through 
the Universe interact with CMB photons producing pi- 
ons, if their energy is above the single pion production 
threshold, 

ml + 2m„M s . , 

^th = -: ■ (5) 

4e 

Here, and Ms denotes the mass of a pion and of 
S, respectively. For the following simple estimates, we 
neglect the Bose-Einstein distribution of the photon en- 
ergies and use just the average energy of a CMB photon, 
e = 7r 4 T /30C(3) w 6.4 x 10~ 4 eV. In order to avoid 
the GZK cutoff, the cross section of strongly interacting 
hadrons S with photons should be smaller than the one 
of nucleons, 

cr S 1 = B \ ° V1 ' 

where the suppression factor in the resonant case comes 
from the different center mass momenta in the Breit- 
Wigncr formula. The dimensionless parameter B is 
B(s) < 1, because the assumed resonant state (the equiv- 
alent of the A resonance) has a mass larger than S. 

Thus, strongly interacting UHE particles with E 3> 
E t h will interact with CMB photons on the typical scale 

hat = _J_ =8Mpc^l, (7) 
a Sl n a Sl 
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where the CMB number density is no = 410 cm -3 and 
cr p7 = 10~ 28 cm 2 is the multi-pion production cross sec- 
tion. During each interaction, the particle S loses the 
fraction y ~ 0.5 of its energy until its energy is close to 
E t h- There, the energy fraction lost reduces to m^/Ms, 
while the cross section can be increased due to reso- 
nances. 

If one defines the radius -Rgzk as the distance after 
which a particle S with E 3> E t h lost 95% of its initial 
energy, then 



R 



GZK 



In 20 

ya Sl n, 



48 ^iMpc . 



(8) 



In the case of protons, this effect was first considered by 
Greisen, Zatsepin and Kuzmin in 1966 and is called 
GZK effect. The threshold energy for proton from Eq. (JJJ 
is E = 1 x 10 20 eV. However, protons with energies 
Egzk = 4 x 10 19 eV can still interact with the high- 
energy tail of the CMB distribution. 

Let us now consider the case of a new, strongly inter- 
acting particle S from a general point of view. If its heav- 
ier than a proton, then its GZK cutoff is both softened 
and shifted to higher energies. The first effect arises be- 
cause of the smaller energy transfer near threshold, while 
the second one is due to a smaller resonant cross section 
with CMB photons. Let us now turn to our specific ex- 
amples G = (gg), and B° = (Id) and B ± = (bu). The 
first case was studied in Ref. 43] . It was found that erg, 
is at least a factor 8 smaller than <r P7 , even for as low 
masses as rrig = 1.5 GeV. This small cross section leads 
together with the reduced energy losses per scattering to 
a shift of the GZK cutoff close to the maximal energies 
of astrophysical accelerators . While in the case of G 
some information about the mass spectrum of low-lying 
g containing hadrons is available [85j , this information is 
missing for B Q ± . Since the knowledge of the low- lying 
resonances of B is essential to perform a detailed calcu- 
lation of its energy losses on CMB photons, we can only 
estimate the energy losses. To be conservative, we as- 
sume that the resonant contribution to Bj scattering is 
only suppressed by its larger mass, or B = 1 in Eq. ||§J, 
and by the smaller energy transfer close to threshold, 
y = m n /Mg. Using as smallest value for Mg ~ 2 GeV 
then shifts E t h by a factor two, what together with the 
other suppression factors causes only a mild GZK effect 
at high enough energies. The resulting spectrum is shown 
for B°, an injection spectrum E~ 2 and uniformly dis- 
tributed sources for a Mg = 2 and 3 GeV in Fig. [21 In 
the case of the charged B ± 1 additionally e + e~ pair pro- 
duction has to be considered. For comparison we show 
in Fig.[3]the proton spectrum from the same distribution 
of sources. 

Another important condition is that the particle S 
should be stable, traveling through the Universe. The 
lifetime ts should be bigger then 
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FIG. 3: Energy spectrum of B° hadrons with injection spec- 
trum E~ 2 and uniformly distributed sources for a Mg = 2 
and 3 GeV; for a comparison proton spectrum is also shown. 
At energies 10 18 eV < E < 6 x 10 19 eV proton spectrum also 
suppressed due to e + e~ production. 



where R is measured in Gpc, Ms in GeV and E20 = 
E/(10 w eV). Note that this allows the possibility that 
the gravitino is the LSP and the gluino or bottom squark 
is the NLSP. In this case, the NLSP decays only via 
gravitational interaction and thus has a long enough life- 
time to serve as UHECR messenger. On the other hand, 
all experimental constraints from searches for anomalous 
heavy isotopes can be avoided easily in this scenario. 

Thus any new strongly interacting "messenger" parti- 
cle with multi-GeV mass and lifetime bigger than a year 
can travel over cosmological distance and solve the GZK 
problem. In particular, gluino and bottom squark con- 
taining mesons and baryons can serve as as "messenger" 
particles. 



V. INTERACTIONS IN THE ATMOSPHERE 

The interactions of glueballinos with nucleons were 
considered in detail in Ref. There, the Monte Carlo 
simulation QGS JET |8(| , which describes hadron-hadron 
interactions using the quark-gluon string model of the 
supercritical Pomeron in the framework of the Gribov- 
Regge approach, was extended to include light gluinos. 
Moreover, extensive air showers (EAS) were simulated 
and the resulting lateral and longitudinal shower profiles 
were compared to those of EAS initiated by protons. The 
authors of Ref. concluded that glueballinos with mass 
> 5 GeV resemble a penetrating particle and can be al- 
ready excluded using existing data, while EAS initiated 
by glueballinos with mass < 3 GeV can be only differen- 
tiated from proton showers by future experiments with 
larger statistics. 

The calculations of Ref. 01 were done only for the 
special case of a glueballino. However, B-hadrons with 
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the same mass should have very similar interaction prop- 
erties. Main reason for this is that the coupling of the 
Pomeron to a hadron as well as the slope of its coupling 
depends essentially on the size of the hadron, and there- 
fore on its reduced mass. Minor differences arise because 
of the different constituent masses of quarks and gluons 
resulting in different momentum distributions of gluino 
and squarks in different hadrons with the same mass. 
Otherwise, the soft and semi-hard interactions have the 
same dependence on its mass. Finally, the hard inter- 
actions of the constituents at UHE energies are practi- 
cally mass independent in the low mass range of interest. 
We conclude therefore that also i?-hadrons with mass 
< 3 GcV produce EAS consistent with present observa- 
tions. 



VI. PRODUCTION OF SUSY PARTICLES IN 
ASTROPHYSICAL SOURCES 

Protons are the most natural candidates for the ob- 
served UHECR with energies above the ankle E > 
10 19 eV. There are several mechanisms which could be 
responsible for the acceleration of protons to the highest 
energies: The most popular one is particle acceleration in 
shock fronts or Fermi acceleration of the first kind. How- 
ever, there are other, more exotic mechanisms as, e.g., 
particle acceleration in the vicinity of black holes rotat- 
ing in an external magnetic field (see for example |87|'). 

Independent of the specific acceleration mechanism, a 
simple estimate of the maximal possible energy up to 
which a source can accelerate particles was suggested by 
Hillas (88|. It is based on the relation i? max = qBL, 
where q is the charge of the accelerated particle, B the 
magnetic field strength in the acceleration region and L 
its size. Only few astrophysical objects are able to accel- 
erate particles to UHE according to this simple criterion. 
Plausible candidates for acceleration to UHE are AGNs 
and several AGN subclasses were suggested as sources 
of UHECRs. The general perception is that it is pos- 
sible to accelerate protons in objects like AGNs up to 
-E-max < 10 21 eV, but that acceleration to higher energies 
is extremely difficult because of energy losses. 



A. Proton-proton interactions 

We start with the perturbative calculation of the pro- 
duction cross section of bottom squarks in proton-proton 
collisions. The main contribution to the total cross- 
section is given by the gluon^gluon subprocess, first cal- 
culated at leading order in [89|. The center-of-mass en- 
ergy y/s of this process is rather high, y^i > 300 TeV, 
and we restrict ourselves therefore to a leading-order cal- 
culation. We have used the Cteq6 parton distribution 
functions (pdf ) [9(J with the scale fj 2 = s and calculated 
as main contribution to the total production cross section 

the parton subprocesses gg — > gg and gg — > bb. For the 



relatively large gluino masses still allowed, nig > 6 GeV, 
the cross section of gg — > gg is too small and we do not 
discuss this case in the following. 
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FIG. 4: (a) Bottom squarks production cross section in 
proton-proton interaction as function of UHE proton energy, 
(b) Energy fraction transfered to bottom squarks from initial 
protons as function of UHE proton energy 

The sbottom production cross section as function of 
the UHE proton energy E p in the lab frame is shown 
in Fig. for masses Mg = 3,4 and 5 GeV. At ener- 
gies E p ~ I0 13 GeV, this cross section reaches several 
mbarn. The fast increase of the cross section with energy 
is caused by the growing number of accessible soft gluons 
with smaller and smaller x m in(s) values. Therefore, the 
prize to be paid for such a large cross section is the small 
energy fraction transferred, (yi s ) — (E^/E p ) ~ I0 -3 
for E p = 10 20 eV, see Fig. \$Jp. Such small values of 
y make it impossible to produce UHE sbottoms in pp- 
collisions: even if the primary protons would have energy 
E p = I0 23 eV, the average energy of the produced sbot- 
tom would be only I0 19 eV. Since 10 4 more energy will be 
dumped into neutrinos and photons than into sbottoms, 
it is impossible to explain the UHECR flux Fcr, 



F C r(E) = 



I0 2U eV 
E 



eV 



(10) 



with sbottoms without overproducing photons and neu- 
trinos. The produced photons will cascade down to GeV 
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energies and overshot the diffuse gamma-ray flux mea- 
sured by EGRET [34| by two orders of magnitude. In 
principle, one can argue, that it is possible to transfer this 
energy already in the source to energies below those mea- 
sured by EGRET energies, thereby avoiding the EGRET 
bound. However, at the same time the neutrino flux of 
the order of 10 5 eV/(cm 2 s sr) will overshut existing lim- 
its on neutrino flux, given by Fly's Eye [3{| , AGAS A p3(| 
and RICE |33. 



B. Proton-photon interactions 
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FIG. 5: (a) Perturbative part of bottom squarks production 
cross section in proton-photon interaction as function of cen- 
ter mass energy, (b) Part of the energy transfered to bottom 
squarks from initial protons as function of center mass energy 
in perturbation theory. 

We consider next the perturbative calculation of the 
production cross section of bottom squarks in proton- 
photon collisions. Now, the most important subprocesses 

for squark production is the parton subprocess jg — > bb. 
The cross section for sbottom production is shown in 
Fig. G3i. Even compared to the relatively small total jry 
cross section which is of the order of 0.1 mb this cross 
section is small. However, now the center-of-mass en- 
ergy can be much smaller: the typical energy of an in- 
frared/optical target photon is in the range O.f-fO eV, 



hence y/s is between the production threshold and sev- 
eral f 00 GeV. Therefore, the energy fraction transferred 
is now much higher, see Fig. \^jp. However, the combina- 
tion of these two suppression factors is again very small. 
Normalizing the UHE j or t flux to the UHECR flux, 
Eq. (|f 0(1 . would produce 10 6 times higher neutrino and 
photon fluxes, which is again in contradiction both with 
diffuse gamma-ray flux measurements and with neutrino 
bounds. 

Apart from the perturbative contribution to the cross 
section calculated above, non-perturbative contributions 
have to be considered where at small momentum transfer 
hadrons interact with each other. A calculation of this 
kind could be performed in the vector dominance model. 
Then the total production cross section can be split into 
two parts, 

a pl » ^ a Wi a w + — a^ qq (,s, fej 

i "i^^o b 

(11) 

where the sum i extends over the vector mesons i with 
weight Wi and the second part describes the perturba- 
tive process 7 — > bb matched to the first contribution at 
Q 2 = Qo HH- The second contribution can be evaluated 
at UHE and gives a pl ~ (rn^/m^apj. The dominant 
subprocess of the first part is the t channel exchange of 
B mesons. It is therefore natural to expect that this con- 
tribution is also suppressed relative to the photon-proton 
total cross-section in the multi-pion production region by 
the ratio m 2 /Mj. We shall parameterize therefore the 
bottom squark production cross section as 



where A(s) < I is a dimensionless factor depending on 
s. We expect A(s) ~ 1 in the high-energy region and 
A(s) — ► for s — >■ 4m?. The transferred energy can be 
as high as f — 50%. The required photon energies are of 
the order of O.f - fO eV. 

Since we are interested in shadron masses around f .5- 
3 GeV, we can compare the value predicted by Eq. i|f 2|) 
with the total charm production cross section in proton- 
photon collisions. Reference [92J collected experimental 
data 93] in the energy range from \/s = 1 0-200 GeV; 
the cross section increases from f0~ 3 mbarn at y/s = 
20 GeV to 10" 2 mbarn at y/s = 200 GeV. According to 
p2jl. we would expect a cross section (to w /Md) 2 ct P7 ~ 
5 x i0~ 3 (7 P7 ~ f0~ 3 mbarn. Hence we conclude that 
Eq. I|f 2|) is a rather conservative estimate. 

Close to threshold, t channel exchange of B mesons 

proceed as p + 7 — > (udb) + (bu) and p + 7 — > (uub) + (ud) . 
Thus at moderate UHE energies, the UHECR flux should 
consist of the usual protons, and positive and neutral b- 
hadrons. At the highest energies, when several 6-hadrons 
are produced, additionally negatively charged 6-hadrons 
appear. 
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Let us compare these numbers with the parameters 
of astrophysical objects. The important difference to 
proton-proton interactions is that the energy of back- 
ground photons and as a result the center-of-mass energy 
is normally much smaller than the proton mass. The re- 
quired center-of-mass energy to produce particles with 
mass in the multi-GeV range is around s = 50 GeV 2 sso. 
It should be somewhat higher than the threshold 4M 2 to 
avoid the kinematical suppression effects near threshold. 
The typical photon energy then is 



2E„ 



= 0.25^- eV. 



3 



20 



(13) 



Photons of such energies exist in many astrophysical 
objects which can accelerate protons. However, the accel- 
erated protons should interact inside these objects with 
photons. In other words, the propagation length i; nt of 
protons should be smaller than the size R of the interac- 
tion region, 



^int — 



1 



< R. 



(14) 



For example, the time variability of several days in the 
optical spectrum of AGN cores corresponds to a region 
size of R = 10 16 cm [9^]. The photon density can be 
estimated from the optical/infrared luminosity, 
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cm 6 R( 6 e-i 



(15) 



the parameters of the astrophysical object are the same 
as in Eq. I|l(j|) . In the case of a light bottom squark we 
suppose that 0.1 is also a reasonable estimate for the 
suppression factor in the <7s 7 cross section and leave the 
detailed analysis for future investigation. 

Thus, new light hadrons can be produced in astrophys- 
ical objects from 10 21 eV protons, interacting with in- 
frared/optical photons of the energies 0.1 — 10 eV, if the 
sources are optically thick for protons, Eq. I|16l) . The 
model for such a source can be similar, for example, to 
the one of Stecker et al. [9(|. Produced hadrons will 
escape from the same objects if their interactions with 
photons are suppressed in comparison to those of pro- 
tons, Eq. (|17|l . However, simultaneously with the new 
hadrons large fluxes of neutrinos and photons will be 
produced unavoidably. In next section we will discuss 
the experimental constraints on these fluxes. 

We have used AGN cores as a working example of as- 
trophysical accelerators, which, as we have shown, can 
obey the condition of high optical depth for protons, 
Eq. I|16(l . and allow shadrons to escape, Eq. I|17(l . Any 
other astrophysical object, which is able to accelerate 
protons to the highest energies and obeys these condi- 
tions can be a source of shadrons as well. 



VII. ALLOWED PARAMETERS OF NEW 
HADRONS CONSISTENT WITH GAMMA-RAY 
AND NEUTRINO BOUNDS 



where the quantities introduced are the dimensionless lu- 
minosity L44 = i/(10 44 erg/s), the region size i?i 6 = 
i?/(10 16 cm) and the typical photon energy e_i = 
e/(0.1 eV). Substituting the multi-pion production cross 
section <j pl « 10 _28 cm 2 and Eq. (|15fl into the condition 
(|14f> we obtain 



R_ 

km 
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L 4 



i?16£-] 



> 1 



(16) 



Thus, if the parameters of the source are similar to 
those in Eq. Ijltjfl . protons produced inside such sources 
will interact with background photons and can poten- 
tially produce secondary hadrons S. However, the pro- 
duced hadrons still need to escape from the astrophysical 
source. The escape condition is inverse to the one given 
in Eq. I|14|l : The optical depth for the new hadrons should 
be small assuming the same source parameters. Then the 
escape condition is 



1 



(17) 



where the optical depth for protons t p1 is defined by 
Eq. ijTtjl) . In the case of a light glueballino, the sup- 
pression of the as-y cross section can be of the order of 
0.1 in comparison to a PJ at center mass energies of the 
order of 10 GeV H3. This is consistent with Eq. (JT7|>, if 
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FIG. 6: Flux of new hadrons S (thick solid line) and protons 
(dashed line) together with cosmic ray data from AGASA 
and HiRes 6']. Protons accelerated to the energy E = 10 21 
eV (line Pmi) produce secondary photons (dashed-dotted line) 
and neutrinos (dotted line). Photon flux constraint from 
EGRET [3^| and upper limits on the diffuse neutrino fluxes 
from AGASA HH, the Fly's Eye jp, the RICE [37|, and the 
Goldstone experiment (GLUE) [38] as indicated. 



As it was shown in the previous section, only interac- 
tions of UHE protons with infrared/optical background 
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FIG. 7: The neutrino flux for one flavor in the model used 
in Fig. |S] and sensitivities of the currently being constructed 
Auger project to electron/muon and tau- neutrinos |9g|. and 
the planned projects telescope array (TA) |99l| (dashed-dotted 
line), MOUNT 0, and, indicated by squares, OWL [l0(ij ]. 
NT200+Eo3, ANTARES [Ht, AMANDA-II and ICE- 
CUBE 54], as indicated. Also shown (dashed line) is an 
extreme scenario with initial proton spectrum l/E 2 , leading 
to a neutrino flux extending to relatively low energies where 
Baikal, ANTARES and AMANDA-II will be sensitive, and 
the atmospheric neutrino flux for comparison. 



photons of the energy 0.1 — 10 eV can produce a sig- 
nificant amount of new strongly interacting hadrons S, 
without overproducing the diffuse photon and neutrino 
backgrounds. The essential condition for this mechanism 
is the high optical depth for protons, Eq. I|16f) . However, 
this condition has not to be fulfilled by all UHE proton 
sources. Sources with small optical depth for protons 
will just emit UHE protons, which will be responsible for 
UHECR below the GZK cutoff. The few sources of UHE 
shadrons will be only responsible for the highest energy 
cosmic rays with E > 10 20 eV and may be, partly, for 
the clustered component at lower energies. It is also not 
excluded that few sources of this kind will responsible for 
most, or all UHECR above 4 x 10 19 eV. 

We do not specify the proton acceleration mechanism 
in the astrophysical objects here. We just suppose that 
these sources can accelerate protons up to 10 21 eV. The 
proton spectrum should be relatively hard in this case, 

1 / E a with a < 2. The exact shape of the spectrum is not 
important, because the optical depth for protons is high 
and most of them are absorbed. The pion production 
threshold for an optical background of 10 eV is E t h — 

2 x 10 16 eV. Thus the flux of protons with energy above 
this threshold will be reduced by the factor exp(— t(E)). 
If the optical depth is high, t(E) ^S> 1, all protons will 
be absorbed in the source, and the proton flux does not 
overshoot the measured spectrum of UHECRs. 

We have used the propagation code ^63] for the calcu- 
lation of the energy spectra of secondary protons, pho- 
tons and neutrinos. As initial spectrum, we have chosen 



a proton spectrum peaked at the energy 10 21 eV, see 
Fig. El The continuation of this spectrum to lower en- 
ergies is possible for any power law up to a < 2. Even 
an initial proton spectrum with l/E 2 will not contra- 
dict UHECR observations at high energies, but only will 
lead to a higher flux of UHE neutrinos at the energies 
10 16 eV < E < 10 20 eV. 

In Fig. El we present one example of such a calcula- 
tion. Cosmic ray data from AGASA [j| and HiRes [(| 
are shown with errorbars. The contribution of the new 
hadrons S to the UHECR spectrum at the highest ener- 
gies is shown with a thick solid line. We have used conser- 
vatively the AGASA spectrum at the highest energies as 
normalization - choosing HiRes data as reference would 
increase the parameter space for the new hadron. The ex- 
act shape of the S particle spectrum is not well defined, 
because it depends on many unknown parameters like 
the spectrum of background photons in the source, the 
distribution of the sources, initial proton spectrum and 
the energy dependence of the production cross section. 
However, the amplitude of this flux is related to the am- 
plitude of the initial proton flux through Eq. (|12|) . This 
fixes for any given mass Ms, parameter A in Eq. i|12|) 
and average transfer energy (y) = (Es/E p ) the normal- 
ization of the initial proton flux, which is marked by pi n i 
in Fig. El The value of the initial proton flux shown in 
Fig. El corresponds for M s = 2 GeV to (E s /E p ) « 0.1 
and A — 1. 

As background photons in the source we have used in- 
frared/optical photons with energies O.leV < e < 10 eV 
and number density n 7 = 5x 10 12 /cm 3 . This corresponds 
to the luminosity L = 5 x 10 43 erg/s, if the radius of the 
emission region is R = 10 16 cm. After several interac- 
tions with the background photons the accelerated pro- 
tons lose all their energy and produce photons and neu- 
trinos. The neutrino flux should obey the existing exper- 
imental limits of AGASA [H, Fly's Eye H, RICE H3 
and GLUE [3^], which are shown in Fig. El Photons 
cascade down to the GeV and sub-GeV region. The ex- 
isting diffuse gamma-ray flux measurement by EGRET 
restricts the photon flux in the MeV-GeV region, see 
Fig. El However, if high magnetic fields exists in the 
source, then part of the photon energy can cascade down 
into the sub-MeV region, where the bounds on the diffuse 
photon background are at least a factor 10 times weaker 
than in the GeV region. Another part of the photon 
flux can produce large scale jets [9?J, thereby again re- 
distributing energy into the sub-MeV region. This uncer- 
tainty of the value of the photon flux makes the existing 
bounds on the neutrino flux much more important. 

Protons escaping from the source at high energies will 
cascade down to energies below the GZK cutoff and can 
contribute to the observed UHECR spectrum, as shown 
in Fig. El We have supposed that there are no UHECR 
sources within Rqzk — 100 Mpc around the Earth. As 
a result, the proton spectrum has a sharp cutoff below 
10 20 eV, see Fig. EJ Thus, if no nearby UHECR sources 
exist, then even the HiRes data are inconsistent with the 
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minimal model of protons coming from uniformly (but 
rare) distributed UHECR sources. 

In Fig. \7\ we show the UHE neutrino flux (per flavor) 
in our model for two extreme initial proton fluxes: pro- 
tons with a spectrum peaked at E ~ 10 21 eV and for a 
1/E 2 spectrum up to E max ~ 10 21 eV. Both cases are 
consistent with present experimental limits. In the same 
figure, we show the sensitivities of future experiments 
to neutrino fluxes: the Auger project to electron/muon 
and tau- neutrinos 98], the telescope array (TA) [99| . 
the fluorescence/ Cerenkov detector MOUNT 1501 . and, 
indicated by squares, the space based OWL |100| (we 
take the latter as representative also for EUSO), the 
water-based Baikal (NT200+) ANTARES [Eg (the 
NESTOR H3 sensitivity would be similar to ANTARES 
according to Ref. flQ2^ . the ice-based AMANDA-II with 
sensitivity similar to ANTARES and km 3 ICECUBE 0. 
We assume that the proposed water based km 3 detectors 
like GVD [55| and NEMO ^ will have sensitivities sim- 
ilar to the one of ICECUBE. As one can see in Fig. [7| 
future experiments will easily detect UHE neutrino flux 
in any model of new light hadrons. In these models, high 
neutrino fluxes as shown are unavoidable high — contrary 
to the case of neutrino produced by UHECR protons in- 
teracting with CMB photons. In the latter case, the neu- 
trino flux can be as high as in Fig. [7| but could be also 
much lower, depending on the init ial proton spectrum 
and the distribution of sources [l03]. 



VIII. DISCUSSION: BL LACS AS UHECR 
SOURCES. 

The results of the previous section do not depend on 
the particular type of astrophysical accelerator. How- 
ever, we have normalized the astrophysical parameters 
to those of AGN cores by purpose. First of all, AGNs 
are one of the best candidates for proton acceleration 
to UHE. Second, a statistically significant correlation of 
UHECRs with BL Lacs, a subclass of blazars (AGNs with 
jets dire cted to us) with weak emission lines, was found 
|l7l Il04j . Motivated by this correlation, we discuss the 
particular case of BL Lacertae as sources of new hadrons 
in this section. 

In Ref. 0] i it was shown that the correlation with BL 
Lacs requires a new, neutral component in the UHECR 
spectrum. Here we have suggested that this component 
is due to new, neutral shadrons. The trajectories of these 
shadrons should point towards their sources, apart from 
small deflections due to possible magnetic moments. As 
we have showed in the previous sections, shadrons are 
a good candidate for UHECRs and can be produced in 
AGN cores if the condition Eq. (|16[) is fulfilled. Now we 
address the question if th e BL Lacs shown to correlate 
with UHECRs in [H[l0j] obey this condition. 

As an example, we have checked this condition for 
BL Lac RX J10586+5628, which is located at redshift 
z = 0.144 and correlated with AGASA doublet E = 



(7.76, 5.35) x 10 19 eV. First of all, let us note that protons 
with these energies can not reach us from the distance 
D = 680 Mpc (we supposed the actual 'best fit' cosmolog- 
ical model with fl M = 0.3, tt A = 0.7, H = 70km/sMpc). 
The optical magnitude in the V-band of this object is 
15.8, which gives as optical luminosity L = 6 x 10 44 erg/s. 
Since the spectra of BL Lacs are broad in the opti- 
cal region (see, e.g., the spectrum of RX J10586+5628 
in |l05l ]L the density of 0.5 eV photons is high enough in a 
region of size R = 10 16 cm to obey the condition Eq. 116|) 
and thus to produce shadrons from accelerated protons. 

If the UHECR primaries are new particles created in 
proton interactions in the source, large secondary neu- 
trinos and photon fluxes are unavoidable. The neutrino 
fluxes are too small to be detectable by current experi- 
ments, but photons can cascade down into the MeV-GeV 
region in the source, and can be measured. Let us com- 
pare the UHECR flux of the BL Lac RX J10586+5628 
with its gamma-ray flux in the MeV-GeV region, mea- 
sured by EGRET. The two events observed by AGASA 
with energy E = (7.76,5.35) x 10 19 eV allow us to 
estimate the integrated UHECR flux, / dEEF(E) ~ 

0.05eV/cm 2 /s, while the integrated EGRET flux in the 
region 100-800 MeV is approximately / dEEF(E) ~ 
10eV/cm 2 /s [l06|. If we suppose that the EGRET flux 
is mostly due to proton energy losses, the ratio of fluxes 
is 5 x 10~ 3 - a value consistent with our model. How- 
ever, the comparison above can be considered only as an 
order of magnitude estimate: first, the flux measured 
by EGRET could be produced by other interactions. 
Second, the energy injected by protons into electromag- 
netic cascades in the core of RX J10586+5628 can be 
redistributed out of the line-of-sight and thus not con- 
tribute to the EGRET measurement. The next gener- 
ation of TeV gamma-ray tele scop es, like H.E.S.S. |107| . 
MAGIC [HI and VERITAS (Tog, will have sensitivities 
in the 10-100 GeV energy region, allowing to measure 
gamma-ray fluxes from distant sources similar to BL Lac 
RX J10586+5628. Such measurements can be compli- 
mentary to the observations of UHECR from the same 
objects and will allow to restrict or to confirm a wide 
class of UHECR models (including the one we consider- 
ing here) which imply the production of secondary par- 
ticles from protons. 

The high optical depth of photons in Eq. I|16|) guar- 
antees that protons lose energy in the interaction region 
and produce shadrons with a ratio of os/<J pl ~ 5x 10~ 3 . 
For an optical depth of r = 5, only the fraction e~ 5 = 
6.7 x 10~ 3 of initial protons will escape from the source 
without interaction. Thus the flux of produced shadrons 
is in this case similar to the flux of escaping UHE pro- 
tons. This example shows how the same source can be 
a source of UHE protons and at same time a source of 
new UHE hadrons. If the optical depth r is smaller, the 
source will dominantly produce protons, if it is higher, it 
will mostly produce S'-hadrons. 

It will be interesting to check the UHECR data at 
lower energy, E ~ (2 — 4) x 10 19 eV, for correlations 
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with RX J10586+5628. This would be the typical en- 
ergy of protons from this object with z = 0.144 taking 
into account energy losses on the way to the Earth. The 
comparison of the low energy proton flux with a possible 
UHE flux of new hadrons can be used to check the con- 
sistency of our model with the assumption that BL Lac 
are UHECR sources. 

In Ref. |110| . Tinyakov and Tkachev examined corre- 
lations of BL Lacs with the arrival directions of UHECR 
allowing for charges Q = — 1,0, +1 of the primaries. 
They showed that the deflection of charged particles in 
the galactic magnetic field can significantly increase the 
correlation with BL Lacs. If primaries can have charge 
Q = 0, +1, they found that 19 from 57 AGASA events 
correlate with BL Lacs, which have magnitude m < 18 in 
optics. The probability that this correlation is by chance 
is 2 x 10~ 4 . 

They assumed that the charged particles are protons 
and the neutral ones photons. This interpretation has 
two important drawbacks: First, both the highest energy 
event with E = 2 x 10 20 eV and charge Q — +1, and the 
event 16 of Table 1 in (Tl^ with energy E = 4.39 x 
10 19 eV and charge Q = +1, which correlates with a BL 
Lac at z = 0.212, can be only explained as background 
event. Second, they were forced to assume that most of 
the BL Lacs with unknown redshift are located nearby, 

2<0.1. 

Now, if we assume that the UHECR primaries corre- 
lated with BL Lacs are new light hadrons which can have 
charge Q = —1, 0, +1, for example B~ , B° and B + , then 
the assumptions above are not required. Shadrons with 
Q = +1 can easily come from high redshift sources up 
to z ~ 0.5 or even higher. Thus one does not need to 
assume that the BL Lacs with unknown redshift arc lo- 
cated near by n or exclude" unsuitable" sources from the 
Table 1 in |110| . Note also that part of the events with 
Q = +1 can still be protons. Also let us remind that the 
deflection in the magnetic field in the ultra-relativistic 
case does not depends on particle mass and hadrons with 
M = 2-3 GeV with charge Q = +1 will be deflected in 
the same way as protons. 

Moreover, the model with light charged hadrons pre- 
dicts also the existence of particles with negative charge. 
However, because these particles will be produced in the 
sources due to p + j reactions, new particles with Q = +1 
or Q = will dominate. Particles with Q = — 1 can be 
produced in such reactions only well above the produc- 
tion threshold and their expected number should be less 
than the number of particles with Q = +1 and Q = 0. 
Moreover, UHE protons can increase the number of par- 
ticles with Q = +1 at lower energies. Thus, a prediction 
of our model is the existence of a small number of neg- 
atively charged UHE cosmic rays, which average energy 
is larger than events with Q = +1 and Q = 0. 

It is impossible to check statistically the last state- 
ment wi th cu rrent data, however some hint can be found 
in Ref. |104| . The authors of this paper chose as sub- 
set of BL Lacs those which are simultaneously EGRET 



sources. They found that 14 BL Lacs correlate with 65 
UHECR from AGASA and Yakutsk data, if they allow 
as charges Q = +1 and Q = 0. The chance probability 
of this correlation is 3 x 10 -7 which is more than 5<r us- 
ing Gaussian statistics. In this data set 8 BL Lacs out 
of 14 are UHECR sources and emit 13 UHECRs. If one 
suppose, that UHECR primaries with Q = — 1 also ex- 
ist, two more UHECR will correlate with the same BL 
Lacs. These two events have energy E > 5 x 10 19 , which 
is much larger than the average energy in this data set. 
Three other UHECR primaries from this data set can 
have either Q = — 1 or Q = 0. All of them also have 
large energies E > 5 x 10 19 . 

One more assumption made by Tinyakov and Tkachev 
is a cut on t he BL Lac magnitude in the optical range, 
m < 18 [H Ill0| . They found that such a cut maxi- 
mizes the correlation with BL Lacs. However, they were 
not able to explain why the correlated BL Lacs are those 
which are most bright in the optical range. In our model 
of new particle production, such a criterion is obvious: 
the optical background is high enough only in the most 
brightest BL Lacs. Hence, only they are able to pro- 
duce shadrons in p7 reactions. BL Lacs with lower op- 
tical luminosity produce protons, which lose energy and 
contribute to the UHECR spectrum at lower energies. 
Another interesting hint is the value m = 18. In Fig. 3 
of Ref. jfS^] , the dependence of the source magnitude as 
function of redshift was shown under the condition that 
sources are optically thick for protons. This line crosses 
the value m — 18 at redshifts of order z ~ 0.5 — 0.6. 
This distance is similar to the one which shadrons with 
A/ = 2 — 3 GeV can still can propagate. 

Thus, we conclude that the correlation of UHECR with 
BL Lacertae ob je cts which was found in an d investi- 
gated in detail in |l04lll0j suggests that at least some, if 
not most UHECR primaries with £>4x 10 19 should be 
new particles with Q = —1, 0, +1. Explanation of the BL 
Lacs correlation with Q = +1 particles by protons seems 
unlikely The model of new light hadrons, for example 
B~ , B° and B + , naturally explains such correlation as 
well as the cut on the BL Lacs magnitude, m < 18. 



IX. CONCLUSIONS 

The HiRes experiment published recently their 
UHECR data which show a cutoff at the highest ener- 
gies as expected in the conservative model of an uniform, 
continuous distribution of astrophysical sources acceler- 
ating protons up to energies E < 10 21 eV. On the other 
side, a clustered component in the arrival directions of 
UHECR with E>4x 10 19 eV is present in the AGASA 
data with a statistical significance of 4.6a\ If one assumes 
that this clustered component is due to point-like astro- 
physical sources, the predicted total number of sources of 
UHECR with E > 4xl0 19 eVis of the order of 400-1200. 
In this paper we showed that this number of sources is 
so small that the model of continuously distributed pro- 
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tons sources is a bad approximation at the highest en- 
ergies. The latter approximation requires 10 3-4 times 
more sources than estimated from the clustering data. In 
other words, the closest proton source is located outside 
the GZK sphere R > 50 Mpc and the energy spectrum 
of UHECR has a sharp exponential cutoff at the energy 
E < 10 20 eV, which is inconsistent even with the HiRes 
data, see Fig. ^ Including the AGASA data makes this 
discrepancy even worse. 

Moreover, a statistically significant correlation at the 
level of 4er of the arrival directions of UHECR with BL 
Lac objects was found 01 • The closest BL Lacs with 
known redshift are located at cosmological distance, z = 
0.03, and protons with E > 10 20 eV cannot reach us from 
these sources. Some events at lower energies also can not 
be protons, because the redshift of these sources is too 
high. For example, BL Lac RX J10586+5628 is located 
at z — 0.144 or at the distance 700 Mpc. Protons coming 
from this object can have a maximal energy around 2 — 
4 x 10 19 eV, while the correlated UHECRs have much 
higher energies, E — (7.76, 5.35) x 10 19 eV. 

Our findings above suggest the existence of particles 
that can be produced at distant astrophysical objects like 
BL Lacs, propagate through the Universe without signif- 
icant energy losses and produce air-showers in the Earth 
atmosphere similar to those of protons. 

In this work we have investigated the possibility that 
such particles are new light hadrons. We have showed 
that such hadrons can be produced in astrophysical ob- 
jects in interactions of accelerated protons with a back- 
ground of optical photons, if the size of the interaction re- 
gion is larger than the interaction length of protons. The 
interaction of the new hadrons with background photons 
should be suppressed to allow them to escape from the 
sources without significant energy losses. This fact as 
well as the requirement that the energy losses of the new 
particles propagating in the CMB are suppressed com- 
pared to protons restrict the new hadrons to be heav- 
ier than 1.5-2 GeV. Since the primary protons also pro- 
duce also large neutrino and gamma-ray fluxes, which 
are bounded by experimental limits and measurement, 
only hadrons with masses below 3 GeV are allowed. The 
possibility to travel over cosmological distances without 
decay restrict the lifetime of these particles to be larger 
than one month. 

As specific example we have considered hadrons con- 
taining light bottom squarks. This case agrees with all 
existing astrophysical observations, if the shadron mass 
is in the window 1.5 GeV < Ms < 3GeV. Such a new 
hadron can explain the observation of UHECR at the 
highest energies. 

If BL Lacs are indeed UHECR sources, our model of 
new light hadrons allows to solve several puzzles con- 
nected with these objects. First, all correlated UHECR 
with zero charge can be our new hadrons with Q = 0. 
Second, our model offers a simple explanation why only 
optically bright BL Lacs correlate with UHECRs: only 



if the density of optical photons in a BL Lac is high, 
the probability of protons to interact and to produce our 
new hadrons is large enough. Magnitude m — 18 can 
correspond to the redshift z ~ 0.5 — 0.6, a distance from 
which our new particles can still reach the Earth without 
significant energy losses. 

In Ref. |110| it was shown that the correlation with 
BL Lacs increases if one supposes that some UHECR 
have non-zero charge. In particular, a significant corre- 
lation was found if some UHECRs have a positive charge 
Q = +1. It was suggested that these positively charged 
particles are protons. However, this assumption forced 
the authors of Ref. |110| to assume that most of the BL 
Lacs with unknown redshift are located at the distances 
z < 0.1. Furthermore, they had to assume that some of 
the UHECR which can not be protons are correlated just 
by chance. 

These two assumptions can be relaxed in our model if 
one assumes that new hadrons with non-zero charge are 
also long-lived. However, the existence of such hadrons 
is disfavored from accelerator experiments. 

An important consequence of our model is an unavoid- 
able high UHE neutrino flux. This flux is well within 
the sensitivity region of all future UHECR experiments 
and can also be detected by km 3 neutrino telescopes like 
ICECUBE (or GVD and NEMO). In the case of an ini- 
tial proton spectrum cx 1/E a with a ~ 2 even 0.1 km 3 
neutrino telescopes like AMANDA II, ANTARES and 
NESTOR will be able to detect the diffuse neutrino flux 
in the 10 16 eV energy region. 

Another consequence of our model is a cutoff in the 
UHECR spectrum, which can be observed around E ~ 
10 21 eV at future UHECR experiments like the Pierre 
Auger Observatory, the telescope array and EUSO. 

New hadrons with 1.5-3 GeV mass can be searched 
for in existing accelerator experiments like CLEO and 
B-factories or with a dedicated experiment as proposed 
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